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FOREWORD 


This  report  was  prepared  by  Physics  International  Company, 
2700  Merced  Street,  San  Leandro,  California  94577,  under 
Contract  No.  F086 35-72-C-0229  with  the  Air  Force  Armament 
Laboratory,  Eglin  Air  Force  Base,  Florida.  Mr.  Lovonia  J. 
Theriot  (DLYA)  managed  the  program  for  the  Armament  Laboratory. 
This  effort  was  conducted  during  the  period  from  June  1972  to 
April  1973. 

The  contractor  report  number  assigned  is  PIFR-430. 

This  report  is  divided  into  two  volumes.  Volume  I  presents 
the  generalized  analytical  approach  to  shaped-charge  warhead 
design.  Volume  II  describes  the  modification  and  utilization  of 
a  two-dimensional  finite  difference  continuum  mechanics  code 
utilizing  the  Lagrangian  coordinate  system  to  calculate  the 
complete  jet  formation  parameters  for  any  generalized  axisym- 
metric  shaped  charge.  This  is  Volume  I. 

This  technical  report  has  been  reviewed  and  is  approved. 

Qj-id. 

ROBERT  W.  DILLON,  Colonel,  USAF 
Chief,  Weapons  Systems  Analysis  Division 


ABSTRACT 


This  report  describes  a  technique  to  optimize  the  current 
shaped-charge  design  procedure  as  follows.  Starting  with  the 
desired  target  to  be  defeated,  a  determination  of  the  desired 
penetration  characteristics  of  t.1  e  jet  would  be  made.  Existing 
jet  penetration  theory  would  then  be  used  to  estimate  the  ideal 
characteristics  of  the  jet  to  defeat  the  given  target.  A  shaped- 
charge  launcher  would  then  be  designed  to  give  these  ideal  jet 
characteristics.  However,  a  suitable  design  procedure  requires 
(1)  a  viable  analytical  or  empirical  design  approach  to  obtain  a 
first  cut  shaped  charge  design,  (2)  a  better  understanding  than 
now  exists  of  the  detailed  mechanisms  of  jet  formation,  and  (3)  a 
better  understanding  of  the  phenomenon  of  jet  penetration-  This 
report,  which  is  contained  in  two  volumes,  addresses  the  first 
two  of  these  requirements.  Volume  I  describes  the  use  of  the 
existing  non-steady  state  theory  of  jet  formation  with  experi¬ 
mental  data  and  one-dimensional  finite  difference  continue  mech¬ 
anics  calculations  to  obtain  the  liner  collapse  velocity  for 
generalized  axisymmetric  shaped  charges.  The  results  of  this 
work  are  then  used  to  obtain  nonunique  shaped  charge  designs  which 
give  the  required  idealized  jet  parameters.  Volume  II  describes 
the  modification  and  utilization  of  a  two-dimensional  finite  dif¬ 
ference  continuum  mechanics  code  utilizing  the  Lagrangian  Co¬ 
ordinate  system  to  calculate  the  complete  jet  formation  parameters 
for  any  generalized  axisymmetric  shaped  charge.  The  utilizatioii 
of  this  code  allows  a  more  detailed  study  of  such  phenomena  as 
jet  stability,  bifurcation  on  the  axis,  shear  gradients,  viscosity 
shocks,  incipient  vaporization,  surface  tension,  and  possible 
other  effects.  The  combined  use  of  both  the  engineering  formula¬ 
tions  along  with  the  sophisticated  two-dimensional  code  calcula¬ 
tion  allows  design  engineers  the  versatility  to  design  the  most 
optimum  shaped  charge  for  their  particular  application. 
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SECTION  I 
INTRODUCTION 

Conical  metallic-lined  shaped  charges  have  been  in  existence 
for  30  years  and  during  this  time  the  basic  design  has  changed 
little.  The  qualitative  theory  of  jet  formation  and  jet  penetra¬ 
tion  was  first  published  in  1948  (Reference  1) .  This  theory  was 
based  upon  the  classical  hydrodynamic  theory  of  perfect  fluids 
and  thus  assumed  steady  state  conditions  to  exist  and  required 
material  flow  with  average  properties  which  conserved  mass, 
momentum,  and  energy.  In  order  to  partition  the  energy  unequi- 
vccably,  the  authors  had  to  ignore  the  role  played  by  shocks  in 
the  acceleration  of  the  conical  liner.  The  theory  presented  in 
Reference  1  assumed  that  each  element  of  the  conical  liner  col¬ 
lapsed  with  a  constant  velocity  and  thus  the  resulting  jet  was 
one  with  uniform  properties.  Uniform  properties  imply  that  the 
jet  mass  per-unit  length  of  jet  and  the  jet  velocity  for  each 
element  of  the  jet  were  equal  for  the  total  length  of  the  jet. 

The  steady  state  theory  of  jet  formation  was  modified  in  1952 
(Reference  2)  to  account  for  the  variable  liner  collapse  velocity. 
This  nonsteady  theory  for  jet  formation  was  then  re-examined  in 
1954  by  Eichelberger  (Reference  3)  and  has  seen  no  further  basic 
modifications  since  that  time.  Since  1954  the  majority  of  the 
literature  discussing  jet  parameters  has  been  the  correlation  of 
extensive  experimental  data  with  the  nonsteady  jet  formation 
theory  to  obtain  empirically  the  liner  collapse  velocities.  Once 
the  liner  collapse  velocities  were  available  for  a  given  charge 
design  and  explosive,  then  quick  calculations  could  be  performed 
to  determine  the  jet  parameters  from  other  explosive  charges  or 
scaled  versions  of  the  same  design.  However,  if  the  charge-to- 
mass  ratio  (c/m)  of  the  design  of  the  liner  angle  or  geometry 
were  changed  then  additional  experimental  data  was  necessitated 
to  obtain  the  collapse  velocity  as  a  function  of  liner  element 
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for  the  new  design.  The  missing  link  to  the  straightforward 
-ilization  r  the  jet  formation  theory,  as  given  in  Reference  3, 
is  an  equation  or  data  that  gives  the  collapse  velocity  of  the 
line,  as  a  function  of  the  explosive  charge  and  the  radius  of  the 
liner  from  the  axis  of  the  shaped  charge  configuration.  Analyti- 
'  -  ..xpressions  for  the  collapse  velocity  of  flat  plates  as  a 

..unction  of  the  above  parameters  have  been  investigated  by  several 
authors;  two  of  the  most  recent  are  given  in  References  4  and  5. 
Similar  analytcal  solutions  to  the  problem  of  the  implosion  of 
infinite  cylinders  by  explosive  are  not  possible  at  this  time. 

The  basic  reason  for  this  is  that  the  material  at  the  inside  of 
the  cylindrical  liner  experiences  an  acceleration  due  to  the  con- 
V' .gent  pressures,  and  thus  the  velocities  through  the  liner 
chix'kness  are  highly  nonuniform. 

Existing  experimental  data  on  the  collapse  velocity  of 
shaped  charges  plus  one-dimensional  calculations  uo  simulate  the 
collapse  of  a  shaped-charge  liner  are  used  to  develop  a  predictive 
capability  of  the  collapse  velocities  of  imploding  cylinders  as  a 
function  of  c/m  and  distance  of  the  liner  from  the  impact  axis. 
This  empirical  data,  combined  with  the  nonsteady  state  theory  of 
jet  formation,  provides  a  convenient  cool  to  design  a  shaped 
charge  for  given  jet  parameters.  This  design  would  then  consti¬ 
tute  an  initial  geometry  configuration  which  could  be  fine  tuned 
by  the  use  of  a  two-dimensional  finite  difference  calculation. 


SECTION  II 


GENERALIZATION  OF  THE  NONSTEADY 
THEORY  OF  JET  FORMATION 

I.  DERIVATION  OF  JET  PARAMETERS 

The  theory  of  jet  formation  as  given  in  References  1  through 
3  was  restricted  to  a  liner  with  a  constant  cone  angle.  To  obtain 
the  maximum  flexibility  in  the  design  of  a  shaped  charge  to  give 
desired  jet  parameters,  it  is  desirable  to  allow  a  generalization 
in  the  liner  thickness,  geometric  shape,  and  point  of  initiation. 

From  Reference  2,  the  important  equations  are: 

6 (x)  =  sin"1  [  VQ (x)  cos  e(x)/2  UD]  ;  (1) 

vj  (x)  =  VQ(x)  cos  IctU)  +  <$(x)  -  g  (x) /2]  / sing  (x)/2 

(2) 

dnu(x)/dm(x)  =  sin2  $(x)/2  (3) 

dm  (x) /dm (x)  =  cos2  8(x)/2  (4) 

s 

Referring  to  Figure  1 ,  6  is  the  angle  between  the  direction 
an  element  of  the  liner  travels  after  being  struck  by  the  detona¬ 
tion  wave  and  normal  to  the  liner  surface,  VQ  is  the  velocity  at 
which  the  liner  element  travels  toward  the  axis,  and  a  is  the 
angle  between  the  tangent  to  the  liner  at  a  point  x  and  the  axis 
of  the  liner. 

UD  is  the  detonation  velocity  of  the  explosive;  ,  the 
velocity  of  the  jet  element  formed;  8,  the  angle  between  the 
collapsing  liner  wall  and  the  axis;  e  is  the  angle  between  the 
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Figure  1.  Generalized  Ax i symmetric  Shaped -Cl targe 
Configuration 


explosive  detonation  velocity  vector  and  the  tangent  to  the 

liner;  m  ,  the  mass  of  the  slug;  m. ,  the  mass  of  the  jet,  and 
s  j 

m,  the  mass  of  the  liner.  The  masses  m  ,  m.  and  m  are  each 

s  1 

functions  of  x.  The  mass  m  is  that  part  of  the  mass  of  the 
liner  that  is  included  between  the  top  (apex)  (x  =  0)  of  the 
liner  and  the  plane  perpendicular  to  liner  axis  at  x  =  x.  The 
masses  mg  and  itu  are  the  parts  of  m  that  end  up  in  the  slug  and 
the  jet  respectively. 

As  seen  from  Equations  (1  through  4) ,  once  the  collapse  veloc¬ 
ity  V  and  the  collapse  angle  8  is  determined,  then  the  jet  velocity 

Vj  and  jet  mass  dnu  can  be  obtained.  Following  the  work  of 

Reference  2,  the  collapse  angle  is  obtained  as  a  function  of  the 
liner  geometry  and  collapse  velocity. 

From  Figure  1,  the  following  geometrical  relations  are 
developed: 

tan  (a  -  e)  =  (r  ^  -  D)  /  (x  -  d)  ;  (5) 

T(x)  =  i-  ({x  -  d}2  +  (r  (x)  -  D}2]\*  (6) 

D  * 

r  (x)  =  r^(x)  ~  v0(x)  “  T(x)l  cos  la(x)  +  6  (x)  ]  ; 

(7) 

Z(x)  =  X  +  VQ(X)  [t(x)  -  T(x)l  sin  Otl  (x)  +  6  (x) )  (8) 

r^  is  the  radius  of  the  liner  at  a  point  x;  T  is  the  time 
the  detonation  front  takes  to  reach  a  point  x  on  the  liner;  t  is 
the  time  it  takes  for  a  point  x  on  the  liner  to  reach  a  radius  r, 
and  z  is  the  corresponding  x  coordinate.  The  tangent  of  the 
collapse  angle,  8,  at  any  given  time,  t,  is  by  definition  the 
partial  derivative  of  r  with  respect  to  z  at  constant  t. 
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Therefore,  we  have  from  equations  (7)  and  (8): 


3Z 


3x 

3z 


Sr&  v  t  a  cos  a  _ 

3x  ~  vo  '-,x  ~  3x 


t  cos  A 


T  8_S25_t  +  21  oos 
ax  ax 


R}  +  T  coa  Aj , 


(9) 


3V  SV 

|2.  „  1  +  Vt  Jl^iS-A  +  1°  t  gin  A  .  lo  T  in 

3X  .  o  3X  3X  3X 


(10) 


.  Vq  {T  ^  +  ||  8in  A  j 


where  A  =  a  +  5. 

Using  the  following  relations  anu  rearranging  equations  (9) 
and  (10) ,  we  have: 


~  c.os  A  *  -  sin  A  (o'  +  6'); 
~  sin  A  «  cos  A  (a*  +  6'); 


where  the  prime  indicates  differentiation  with  respect  to  x. 


3r 

32 


3x 

3z 


=  M 

3z 


3x 


-  (t  -  T) 


1  +  (t  -  T) 


fr. 


{v°°3 


A  -  VQ  (o*+6 ' )  sin  A^+  VQ  T’ 


in  a|- 


(11) 


(a '+6')  cos  A  +  V  *  sin  A 

o 


j  -  VQ  T'  si 


sin  A 
(12) 


cos  A 


The  collapse  angle  of  interest  is  at  r  =  0  when  an  element  at  the 
point  x  on  the  liner  reaches  the  axis.  Therefore  from  equa¬ 
tion  (7)  ,  the  corresponding  time  at  r  =  0  is  given  as: 


t  -  T  = 


V  COS  A 
o 


Substituting  (13)  into  (11)  and  (12)  and  dividing  (11)  by  (12) 
gives  tan  B  =  3r/3z|  ; 


r  V  •  /V  +  V  T ' 

tan  B  =  tan  a  ±  *  L  (ot 1  -t<S  * )  tan  A  -  o  '  oj  o  cos  A 

1  +  r,  l  (a'+6')  +  V  '/V  tan  A]  -  V  T*  sin  A 
_ _ *  o  o  o 


Therefore,  once  a  liner  geometry  is  chosen  and  the  liner  collapse 
velocities  determined,  the  jet  parameters  can  be  calculated. 

Before  proceeding  with  the  use  of  these  equations,  the 
equations  to  calculate  6 ' ,  e ' ,  and  T '  will  be  obtained . 

From  equation  (1) , 


sin  6  =  VQ  cos  e/2  Ud; 


6' 

cos  6  = 

o 

rv  ‘ 

6' 

=  tan  6 

o 

V 

D 

tan  e 
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The  derivative  of  e  with  respect  to  x  is  obtained  from 
equation  (5)  and  is  given  by  the  following  equation. 


o’  +  SSS-j-j-g-O.  [tan  (a-e)  -  tan  a] 
x-a 


The  derivative  of  T  with  respect  to  a  is  obtained  from 
equation  (6)  and  is  given  by  the  following  equation. 


x  -  d 


UD  T 


[1  +  tan  (a-e)  tan  aj 


Equations  (1  through  5)  and  (14  through  17)  constitute  a 
sufficient  set  to  calculate  the  jet  parameters  from  a  generalized 
axisymmetric  shaped  charge  once  the  liner  geometry  and  collapae 
velocities  are  known. 

2.  EXAMINATION  OF  THE  JET  PARAMETER  EQUATIONS 

Before  proceeding  to  a  determination  of  the  liner  collapse 
velocity#  it  is  beneficial  to  evaluate  the  jet  parameter  equa¬ 
tions  in  more  detail. 

Unfortunately#  the  interdependence  of  the  liner  geometry 
variables  and  the  liner  collapse  velocities  on  the  collapse 
angle  8  do  not  permit  the  determination  of  a  unique  shaped 
charge  design  for  given  jet  parameters.  As  a  result#  an  initial 
liner  geometry  is  chosen  and  a  number  of  calculational  iter-  ■ 
ations  are  made  to  determine  the  required  final  geometry  and/or 
required  collapse  velocities.  To  reduce  some  of  the  preliminary 


iterations,  the  dependence  of  the  collapse  3  on  the  basic 
shaped  charge  variables  have  been  plotted  in  the  following 
figures . 

The  variation  of  the  ratio  of  jet  velocity  to  liner 

collapse  velocity  versus  3  is  shown  in  Figure  2  as  a  function 

of  a  +  6.  Figure  3  shows  the  variation  of  3  with  V'  /V  as  a 

o  o 

function  of  a  and  r.  for  fixed  values  of  V  =  0.25,  6  ~  10°, 

X*  o 

6 '  ~  0 ,  a '  =  0 ,  and  T '  =  1.5.  For  most  shaped  charges  of 
interest,  5  varies  from  5°  to  15°,  6'  is  small,  a'  is  zero  for 
a  constant  angle  liner,  and  T'  is  approximately  1.5.  Thus,  this 
figure  provides  a  good  approximation  for  a  constant  angle  liner 
geometry  and  can  be  useful  in  selecting  the  liner  angle.  As  the 
liner  radius  increases,  the  collapse  angle  3  becomes  more 
sensitive  to  the  velocity  gradient  along  the  liner.  This  implies 
that  for  a  nonconstant  collapse  velocity,  the  uncertainties  in 
the  jet  parameters  become  greater  with  increasing  liner  radius 
when  a  is  constant.  It  is  also  noted  that  if  the  collapse 
velocity  is  constant  along  the  liner,  then  the  collapse  angle  is 
also  constant.  The  slope  of  the  lines  shown  in  Figure  3  are  equal 
for  all  collapse  velocities  with  only  a  change  in  position  along 
the  v'0/v0  =  0  coordinate.  The  variation  of  3  with  V  and  a  for 
V'0/Vo  =  0  is  shown  in  Figure  4.  The  combined  use  of  Figures  3 
and  4  will  yield  additional  Figure  3's  for  other  collapse 
velocities. 

The  use  of  shaped  charge  liners  with  a  variable  a  provides 
both  more  flexibility  and  complexity  to  the  problem.  Figure  5 
shows  the  variation  of  3  with  a '  for  other  selected  parameters. 

It  is  interesting  to  note  from  this  figure  that  for  a’  >  3,  3  is 
nearly  independent  of  V0'/V0  and  r^.  and  approaches  the  values 
given  in  Figure  4.  It  therefore  appears  that  a  liner  with  a 
variable  a  will  provide  a  more  reproducible  jet  as  well  as  more 
flexibility  in  the  design. 
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Figure  S.  Variation  of  8  with 


SECTION  III 


DETERMINATION  OF  LINER  COLLAPSE  VELOCITIES 

The  previous  analytical  jet  parameter  equations  assumed 
that  the  shaped-charge  liner  could  have  a  variable  collapse 
velocity  along  its  length  but  not  through  its  thickness.  The 
last  part  of  this  assumption  implies  that  the  stresses  within 
the  liner  are  small.  This  is  far  from  true  and  in  fact  both 
the  velocity  and  stress  through  the  liner  thickness  may  vary 
by  factors  of  two  or  more.  Therefore,  the  collapse  velocity 
used  in  the  previous  equations  must  be  viewed  as  an  equivalent 
steady  state  collapse  velocity  which  may  or  may  not  be  corre¬ 
lated  with  the  actual  collapse  velocity  profile  through  the 
liner  thickness.  The  usefulness  of  the  analytical  jet  para¬ 
meter  equations  is  then  dependent  upon  the  availability  of 
steady  state  collapse  velocity  data  for  imploding  cones  or 
cylinders  for  a  variety  of  liner  masses,  internal  radii, 
charge  mass,  etc.  Since  an  extensive  experimental  and  analyti¬ 
cal  program  (References  6  to  9)  was  required  to  obtain  these 
data  for  the  105  mm  unconfined  shaped  charge,  it  is  impractical 
to  obtain  the  required  data  in  the  same  manner. 

The  approach  taken  in  this  report  to  obtain  the  equivalent 
steady  state  collapse  velocity  ia  to  utilize  the  existing  105  mm 
unconfined  shaped  charge  experimental  data  with  one-dimensional 
calculations  of  the  explosive  implosion  of  cylinders. 

1 .  Approach 

Over  the  past  nine  years,  the  contractor  has  been 
utilizing  linear  explosive  drivers  (References  10  to  13) 


to  launch  projectiles  to  hyper velocity,  simulate  nuclear  blast 
waves,  dr_ve  high  performance  shock  tubes,  plus  other  applica¬ 
tions.  The  linear  explosive  driver  is  simply  a  nonjetting 
cylindrical  shaped  charge.  As  a  result  of  the  extensive  and 
variable  use  of  these  drivers,  the  contractor  has  calculated  the 
collapse  process  with  two-dimensional  f i.  ‘ te  difference  continuum 
mechanics  codes,  approximated  the  collapse  with  one-dimensional 
codes,  and  used  flash  X-ray  photography  to  experimentally  deter¬ 
mine  the  collapse  geometry.  On  the  basis  of  this  work,  it  has 
been  found  that  segments  of  a  conical  imploded  system  can  be 
assumed  to  act  as  an  equivalent  infinite  cylinder  in  a  similar 
geometry.  On  the  basis  of  this  assumption,  one-dimens .ional  finite 
difference  calculations  were  performed  which  simulated  different 
segments  of  the  105  mm  unconfined  shaped  charge  for  which  semi- 
empirical  collapse  velocity  data  exists.  From  a  comparison  of 
the  calculated  velocity  history  with  the  data,  it  was  hoped  a 
correlation  could  be  obtained  that  would  allow  a  determination 
of  the  equivalent  steady  state  collapse  velocity  from  additional 
one-dimensional  calculations . 

2 .  Procedure 

The  BRL  105  mm  unconfined  shaped  charge  (Figure  6)  was  used 
as  the  standard  for  correlation  with  the  one-dimensional  cylin¬ 
drical  implosion  calculations.  The  copper  liner  was  then 
approximated  as  independent  cylindrical  segments,  driven  by  a 
volume  burned  composition  B  explosive.  The  volume  burn  code 
option  sets  the  initial  pressure  in  each  HE  zone  equal  to  one- 
half  the  Chapman- Jouquet  pressure  which  has  been  found  to  corre¬ 
late  well  with  experiment.  A  possible  explanation  for  this 
correlation  is  than  in  a  sideward  burn  the  kinetic  energy  in 
the  detonation  products  behind  the  burn  front  does  not  partici¬ 
pate  in  the  implosion  of  the  cylinder  and  thus  only  one-half  of 
the  HE  energy  is  available.  One-dimensional  calculations  were 
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then  run  for  eleven  cylindrical  segments  at  approximately  0.5  cm 
intervals  along  the  axis  ranging  from  4  cm  to  8.85  cm  from  the 
cone  apex.  The  copper  liner  was  divided  into  fifteen  equal 
radial  zones  and  the  velocity  and  kinetic  energy  of  each  zone 
versus  time  were  monitored.  An  example  of  the  velocity  time 
profiles  at  several  axial  locations  are  given  in  Figures  7  to  14. 

A  review  of  Figures  7  to  14  shows  a  number  of  important 
phenomena.  First,  the  acceleration  of  the  liner  is  discontinuous 
in  time  as  well  as  variable  through  the  liner  thickness..  The 
regions  of  constant  velocity,  or  zero  acceleration,  correspond 
to  the  round  trip  transit  time  of  the  shock  between  the  point 
within  the  liner  and  either  the  outside  or  inside  surface  of  the 
liner.  The  effect  of  the  radial  convergence  of  the  liner  is 
seen  as  a  rapid  increase  in  the  velocity  of  the  inside  of  the 
liner  just  prior  to  complete  collapse.  As  the  liner  radius  is 
increased,  the  velocity  gradient  through  one  liner  at  the  time 
of  collapse  is  reduced.  This  results  from  both  the  reduction 
in  the  amount  of  HE  driving  the  liner  with  increasing  radius 
as  well  as  more  time  for  the  stresses  within  the  liner  to 
equilibrate,  it  is  interesting  to  note  that  at  7  usee,  the 
liner  has  between  86  and  94  percent  of  its  total  kinetic  energy. 

Since  the  liner  collapse  process  is  very  nonsteady,  the 
choice  of  an  equivalent  steady  state  velocity  to  use  in  the 
jet  parameter  equations  will  be  somewhat  arbitrary.  The  assump¬ 
tion  is  made  that  the  liner  has  collapsed  at  the  time  the  velo¬ 
city  of  the  second  zone  has  reached  a  maximum.  The  time  of 
collapse  of  the  first  zone  is  not  used  because  the  zone  size  of 
the  inside  zone  will  limit  the  code  accuracy  during  the  final 
convergence  and  collapse  on  the  axis.  Figure  15  shows  the 
velocity  profile  through  the  liner  at  the  time  of  collapse  at 
various  points  along  the  liner.  Also  shown  are  the  semi-empirical 
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Figure  8.  Collapse  Velocity  Versus  Time  at  R  =  1.535  cm 
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Figure  31.  Collapse  Velocity  Versus  Time  at  R  *  2.82  cm 

(Zones  10  through  IS) 
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Figure  12.  Collapse  Velocity  Versus  Time  at  R  *=  2.82  cm 

(Zones  10  through  15) 
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Figure  13.  Collapse  Velocity  Versus  Time  at  R  =  3.21  cm 

(Zones  0  through  3) 
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Figure  14.  Collapse  Velocity  Versus  Time  at  R  *  3.21  cm 

(Zones  12  through  15) 


x*5.85  cm 


x=7. 35  cm 


x*8.35  cm 


NOTES: 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

Mass  thickness,  l  mi/MT 


1.  105  mm  BRL  unconfined  shaped  charge 

2.  Zero  abscissa  is  inside  radius 

3.  I  m.  is  total  mass/length  to  point  i,  MT  is  total 
i  liner  mass/length 

4.  O indicates  data  from  105  mm  shaped  charge  experiments 


Figure  15.  Velocity  Profile  Through  Liner  Thickness 
at  the  Time  of  Collapse 
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steady  state  collapse  velocity  data  points  from  Reference  9. 
Figure  16  shows  the  average  velocity  per  increment  of  the  liner. 
For  example,  the  average  velocity  of  the  first  25  percent  of  the 
liner  is  found  at  the  intersection  of  an  ordinate  through  the 
abscissa  at  0.25  and  the  appropriate  velocity  curve.  The  assump¬ 
tion  is  now  made  that  either  the  actual  or  average  velocity  of 
the  calculated  collapse  velocities  can  be  correlated  with  the 
steady  state  data.  The  form  of  this  assumed  correlation  is  shown 
in  Figure  17.  The  average  velocity  correlation  was  chosen  for 
use  because  of  it  being  less  sensitive  to  changes  in  the  rela¬ 
tive  point  in  the  liner  thickness  and  the  use  of  average  proper¬ 
ties  is  generally  more  meaningful  than  using  the  actual  velocity 
of  some  point  within  the  liner. 


The  curves  shown  in  Figure  17  were  plotted  as  a  function 
of  I^/R*  independent  of  the  liner  mass.  This  implies  that  the 
•location  of  the  equivalent  steady  state  velocity  within  the 
liner  thickness  is  primarily  dependent  upon  the  internal  radius 
of  the  liner  and  the  amount  of  high  explosive  driving  the  liner 
and  not  on  the  liner  mass.  This  assumption  requires  verification. 
Unfortunately,  due  to  the  lack  of  collapse  velocity  r ata ,  the 
only  way  to  verify  this  assumption  is  by  indirect  methods. 
Assuming  the  two-dimensional  finite  difference  code  calculates 
the  jet  parameters  accurately,  then  that  code  can  be  used  to 
evaluate  the  above  assumption.  The  alternate  method  is  to  cal¬ 
culate  the  jet  parameters  of  a  known  shaped  charge  whose  liner 
mass  is  considerably  larger  than  the  105  mm  unconfined  shaped 
charge . 


Using  Figure  17  as  the  key  to  determine  the  .location  within 
the  liner  thickness  where  the  calculated  average  collapse  velo¬ 
city  is  assumed  to  be  the  equivalent  steady  state  collapse 
velocity.  Figure  18  was  constructed  using  the  results  of  35 

* 

Rq  is  the  outside  radius  at  the  HE. 
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Figure  16.  Average  Velocity  Per  Liner  Increment  at  the 
Time  of  Collapse 
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Figure  17.  Proposed  Correlation  Between  Calculated 
and  Semi-Empirical  Steady  State  Collapse  Velocities 
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one-dimensional  calculations.  These  calculations  used  compo¬ 
sition  B  as  the  high  explosive  and  copper  as  the  liner.  To  a 
first  approximation,  other  liner  materials  may  be  substituted 
for  copper  and  the  scaling  procedure  as  given  in  Reference  14 
can  be  used  to  determine  the  effect  of  other  explosive  loads. 

Since  the  calculated  liner  collapse  velocities  are  keyed 
to  the  empirical  data  from  the  105  mm  unconfined  shaped  charge, 
the  data  given  in  Figure  18  is  most  accurate  for  collapse 
velocities  less  than  0.25  cm/ysec  and  for  cone  angles  near 
21  degrees.  However,  for  initial  design  analysis,  these  curves 
should  suffice. 

Figure  18  and  the  jet  parameter  equations  are  sufficient 
to  either  calculate  the  jet  parameters  from  a  given  shaped 
charge  or  to  design  a  shaped  charge  to  give  desired  jet 
parameters . 


SECTION  IV 


APPLICATION  OF  THE  JET  PARAMETER  EQUATIONS 

Referring  to  Figure  18/  the  effect  on  the  collapse  velocity 
from  changing  either  the  HE  thickness,  the  internal  liner  radius, 
or  the  liner  mass  per  unit  length  is  easily  seen.  For  example, 
increasing  the  HE  thickness  such  that  the  ratio  of  R^/Rq  is 
decreased  below  0.4  does  not  increase  the  collapse  velocity 
significantly.  On  the  other  hand,  increasing  the  internal 
radius  of  the  liner  such  that  the  ratio  R^/RQ  is  greater  than 
0.7  will  place  unreasonable  tolerance  limits  on  the  HE  thickness 
to  assure  a  reproducible  collapse  velocity.  Therefore,  in  terms 
of  good  design  practice  one  would  probably  restrict  the  geo¬ 
metrical  range  of  interest  of  R^/Rq  between  the  limits  of  0.4 
to  0.7. 

1.  BRL  105  mm  Precision  Shaped  Charge 

Since  the  finite  difference  calculation  was  performed  fcr 
the  BRL  105  mm  precision  shaped  charge,  the  first  example  of 
the  use  of  the  jet  parameter  equations  tooether  with  Figure  18 
also  used  this  same  shaped  charge.  Figure  19  shows  the  cross 
section  of  this  shaped  charge  and  Reference  15  discusses  the 
details  of  its  fabrication.  The  liner  of  this  shaped  charge 
is  a  remachined  liner  from  a  105  mm  unconfined  shaped  charge. 
Therefore,  the  thinner  liner  of  the  precision  shaped  charge  with 
respect  to  the  105  mm  unconfined  shaped  charge  should  result  in 
higher  liner  collapse  velocities,  higher  jet  velocities,  and 
smaller  jet  mass. 
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The  liner  collapse  velocities  as  a  function  of  axial  distance 
wei^e  interpolated  from  Figure  18  and  plotted  in  Figure  20.  The 
semi-empirical  collapse  velocities  for  the  105  nut:  unconfined 
shaped  charge  are  also  3hown  in  Figure  20  and  as  anticipated  are 
less  than  those  shown  for  the  precision  shaped  charge.  The 
smoothed  values  of  the  collapse  velocity  along  with  its  axial 
derivative  were  then  used  in  Equation  (14)  to  determine  the 
collapse  angle  8.  The  calculated  jet  velocity  and  mass  are 
shown  in  Figure  21  along  with  data  from  the  105  mm  unconfined 
shaped  charge.  As  seen,  the  shape  and  magnitude  of  the  calcu¬ 
lated  curves  are  as  expected  when  compared  with  the  unconfined 
shaped  charge. 

2.  Design  of  a  Near  Constant  Jet  Velocity  Shaped  Charge 

One  possible  application  of  a  special  designed  shaped  charge 
is  where  large  standoff  distances  are  required.  For  example, 
conventional  shaped  charges  cannot  be  used  if  a  25  foot  standoff 
is  required  since  the  jet  would  have  broken  up  into  many  in¬ 
effective  particles  before  reaching  the  target.  Howe^/er,  if  the 
jet  velocities  were  nearly  constant,  then  the  amount  of  jet 
elongation  and  thus  time  of  breakup  could  be  designed  to  accom¬ 
plish  significant  penetration  at  long  standoffs.  If  a  jet  could 
be  made  that  had  a  constant  j*»t  velocity,  then  only  aerodynamic 
effects  would  limit  its  standoff  distance.  However,  since 
fabrication  tolerances  will  not  accommodate  such  a  design,  a 
diverging  jet  velocity  is  required  to  assure  suf  icient  jet 
length  at  the  target  independent  of  nominal  fabrication  tolerances 
For  an  illustrative  example,  consider  an  initial  jet  with  the 
following  properties: 

Initial  Jet  Length  =  4  cm 

Jet  Tip  Velocity  =  0.58  cm/usec 

Jet  Tail  Velocity  =  0.54  cm/usec 


Axial  distance,  cm 


Figure  20.  Liner  Collapse  Velocity  for  the 
BRL  lOSiran  Precision  Shaped  Charge 


Jet  sees,  9*/c* 


The  tip  and  tail  jet  velocities  were  chosen  to  (1)  stay  above  any 
penetration  cut-off  due  to  target  hardness,  (2)  account  for 
fabrication  tolerances,  and  (3)  limit  the  jet  lengthening  to  a 
minimum  and  maximum  of  2  and  25  over  a  25-foot  distance.* 

After  iterating  on  several  liner  geometries,  the  one  shown 
in  Figure  22  satisfies  the  above  requirements.  The  procedure 
used  to  arrive  at  Figure  22  was  to  start  with  an  initial  cone 
angle  and  collapse  velocity  to  give  the  jet  tip  velocity.  The 
cone  angle  and  collapse  velocity  were  then  varied  incrementally 
until  the  desired  jet  velocity  for  the  next  axial  segment  was 
obtained.  A  plot  of  the  resulting  liner  collapse  velocity  is 
shown  in  Figure  23.  Figure  18  was  then  used  to  determine  the 
liner  mass  and  high  explosive  thickness.  This  last  step  allows 
some  degree  of  flexibility  in  the  choice  of  liner  thickness 
variability  and  thus  jet  mass  if  required.  The  final  jet  para¬ 
meters  for  this  example  are  given  in  Figure  24. 


* 

The  designed  jet  lengthening  is  a  factor  of  12.  The  factors 
of  2  and  25  are  for  the  worst  cases  of  the  jet  velocities 
varying  by  ±  3  percent. 
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Figure  22.  Geometry  for  a  Near  Constant  Jet 
Velocity  Shaped  Charge 
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Figure  24.  Jet  Parameters  for  Shaped  Charge 
Geometry  Shown  in  Figure  22 
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SECTION  V 


CONCLUSION 


The  nonsteady  hydrodynamic  theory  of  jet  formation  has  been 
extended  to  include  variable  geometry  liners.  In  addition,  semi- 
empirical  shaped  charge  collapse  data  has  been  combined  with 
one-dimensional  cylindrical  implosion  calculations  to  obtain 
liner  collapse  velocities  as  a  function  of  explosive  thickness, 
radial  distance,  and  liner  mass  (Figure  18).  The  liner  collapse 
velocities  given  in  Figure  18  are  the  equivalent  steady  state 
velocities  required  in  the  theory  of  jet  formation. 

It  is  now  possible  to  calculate  the  jet  parameters  from 
existing  shaped  charges  without  resorting  to  the  use  of  two- 
dimensional  calculations  or  elaborate  experiments.  More  signifi¬ 
cant,  the  combined  use  of  the  collapse  velocity  curves  plus  the 
jet  formation  theory  provides  a  basis  from  which  a  shaped  charge 
can  be  designed  to  give  predetermined  jet  parameters.  Thus 
special  warhead  applications  may  be  practical  or  perhaps  standard 
shaped  charges  can  be  designed  to  perform  the  same  task  more 
efficiently.  Additionally,  the  curves  can  be  ed  to  estimate 
the  variability  in  the  jet  parameters  from  uncertainties  in  the 
fabrication  process. 

Two  examples  of  the  use  of  the  jet  formation  theory  combined 
with  the  collapse  velocity  curves  from  Figure  18  were  given  in 
the  previous  section.  To  verify  the  adequacy  of  this  design 
technique  it  is  suggested  that  a  special  shaped -charge  warhead 
be  specified,  designed,  fabricated  and  tested.  To  the  extent 
that  this  design  technique  can  be  validated,  the  design  of 
future  shaped-charge  warheads  can  be  greatly  simplified. 
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